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Abstract 

The status of our understanding of the mechanisms of color confinement is reviewed, 
in particular the results of numerical simulations on the lattice. 

1 Introduction 

Quarks and gluons are visible at short distances. They have never been ob- 
served as free particles. Search for quarks started in 1963, when first Gcll-Mann 
introduced them as fundamental constituents of hadrons: the signature is their 
fractional charge q = ±1/3, ±2/3. 

The upper limit to the ratio of quark abundance to proton abundance is 
rig/rip < 10~ 27 ^, to be compared with the expectation in the Standard Cos- 
mological Model n q /n p — 10 _12 j2]. The experimental limit on the production 
cross section in nucleon collisions is pQ a < 10~ 40 cm 2 to be compared with the 
expected value a ~ 10~ 25 cto 2 in the absence of confinement. The only natural 
explanation of these small numbers is that these ratios are exacly zero, or that 
confinement is an absolute property due to some symmetry. 

A transition, however, can occur at high temperature to a phase in which 
quarks and gluons are deconfined and form a quark-gluon plasma. 

Big experiments at heavy ions colliders aim to detect such a phase transition, 
even if no clear signature for it is known. 0] A number of theoretical ideas exist 
on the mechanism of confinement, which we shall briefly review below. 

Lattice is a unique tool to investigate the existence of the deconfining phase 
transition, to check the theoretical ideas and possibly to give indications on 
what to look at in experiments. 

2 Lattice investigations 

The partition function of a field theory at temperature T is equal to the euclidean 
Feynman integral extending in time from to 1/T, with periodic boundary 
conditions in time for bosons, antiperiodic for fermions. 
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Finite temperature QCD is simulated on a lattice iV s 3 x N t with N s > N t . 
The temperature T is given by 



T = (2) 

where a is the lattice spacing in physical units , which by renormalization group 
arguments in the weak coupling regime is given by 

o(/9) = ~exp(-/3/&o) (3) 

with (3 — 2N/g 2 ; —bo the lowest order coefficient of the beta function, which 
is negative because of asymptotic freedom. It follows then from eq.(2) that the 
strong coupling region corresponds to low temperatures, weak coupling to high 
temperatures. 

If a deconfining phase transition exists at some temperature, how can it be 
detected, or what is the criterion for confinement? 

For pure gauge theories a reasonable answer exists, which consists in look- 
ing at the static potential acting between a quark and an antiquark at large 
distances: if it is positive and diverging there is confinement by definition. In 
principle this criterion does not insure that no colored particles exist as an 
asymptotic state, but certainly means that heavy quarks are confined. The 
static potential is related to the correlator D(x) of Polyakov lines 

D(x) = (L(x)L(0)) (4) 

as follows 

V(x) = -TIylD(x) (5) 
It can be shown by use of the cluster property that at large distances 

D(x) ~ ccxp(-^|f|)+|(L)| 2 (6) 

\x\ — >oo \ 1 / 

A temperature T c is found in numerical simulations such that 

for T < T c (L) = and hence V(r) = a ■ r (confinement) 

for T > T c (L) ^ 0, V(r) ~ const (deconfinement) 

both for SU(2) and for SU(3) pure gauge theories. (L) is an order parameter 
for confinement, and Z% is the corresponding symmetry. 

Of course no real phase transition can take place on a finite lattice [S] , so 
that the transition from to 1 of (L) is smooth on a finite lattice , and becomes 
steeper and steeper as the size goes large. The steepness is measured by the 
susceptibility xl, 



d 3 x(L(x)LH0) - L(0)tf(0)) (7) 
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which diverges with some critical index 7 at the critical point 

r=(l-|) (8) 

Other relevant critical indices are the index v of the correlation length £ of the 
order parameter 

i (X r~ v (9) 

and the index a of the specific heat 

C v - C° oc r- a (10) 

a, 7 and v identify the universality class and/or the order of the phase transition. 
A weak first order transition is a limiting case a — 1 , 7 =1 and v— 1/d (d the 
number of spacial dimensions i.e. 3). 

The critical indices are determined from the dependence of susceptibilities 
on the spacial size of the system, by use of a technique known as finite size 
scaling||. The result is that for quenched SU(2) the transition is second order 
and belongs to the universality class of the 3d ising model 7 , for SU(3) it is 
weak lrst order [S]- 

In the presence of dynamical quarks Z3 is not a symmetry and therefore 
the Polyakov line cannot be an order parameter. Moreover the string breaks due 
to the instability for production of dynamical quark pairs and the potential at 
large distances is not growing with the distance, even if there is confinement. 
Another symmetry exists at zero quark masses, the chiral symmetry.At T=0 it 
is spontaneously broken , the pseudoscalar bosons being the Goldstone particles 
, but it is restored at T ~ 170 Mev. The corresponding order parameter is 
{ipij>). It is not clear what exacly chiral symmetry has to do with confinement: 
in any case it is explicitely broken by quark masses , and therefore it cannot 
be the symmetry responsible for confinement discussed in sect 1. For a theory 
with Nf = 2, m u — = m, which is a model approximation of reality, the 
situation is schematically represented in fig 1. The critical line T c (m) is defined 
by the maxima of the susceptibilities xl> Xm> Xc„ j which coincide within errors 
[HI E3j , and as an empirical definition the region below the line is assumed to be 
confined, the region above it to be deconfincd. Theoretical ideas are needed to 
understand the symmetry pattern of the system. 

As for the chiral transition a renormalization group analysis and the assump- 
tion that that the pions are the relevant degrees of freedom gives the following 
predictions |llj. If the U(l) axial symmetry is restored below the chiral transi- 
tion, the transition is first order and such is the critical line at to =/= 0. If instead 
the anomaly persists below T c the transition is second order and the critical line 
at to ^ is a crossover. 

3 Theoretical ideas 

A number of theoretical models of deconfincmcnt exist in the literature. 
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Figure 1: The phase diagram of Nf = 2 QCD. 



There is a Gribov model, which is a clever picture of the chiral phase 
transitional- It does not apply, however, to quenched theory. In the spirit 
of the N c — > oo approach the mechanism of confinement should be the same for 
quenched and unquenched. 

Confinement could be produced by the condensation of vortices ^3] The 
model corresponds to a well defined symmetry in 2+1 dimensions and quenched 
theory , but in any case the Z3 symmetry does not survive the introduction of 
dynamical quarks. 

A most appealing idea is dual superconductivity of the vacuum|14|: chro- 
moelectric charges are confined by dual Meissner effect, which sqeezes the chro- 
moelectric field acting between colored particles into Abrikosov flux tubes, in 
the same way as magnetic charges are confined in ordinary superconductors. 
A number of pioneering papers on this mechanism were based on the defini- 
tion and the counting of monopoles. J7| We shall instead concentrate on the 
symmetry patterns involved. Dual superconductivity means that the vacuum 
is a Bogolubov-Valatin superposition of states with different monopole charge 
(monopole condensation) . 

In order to define monopoles a magnetic U(l) gauge symmetry must be 
identified in QCD , which has to be a color singlet if monopoles condense without 
breaking the color symmetry. 

The procedure to identify such magnetic U(l) is known as Abelian Projection 
|15| . We shall present the abelian projection in a form which will prove useful 
for what follows [TTSj. Let G p „ = T l G^„ be the gauge field strength with T l the 
gauge group generators in the fundamental representation, and $ = T l $ z any 
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operator in the adjoint representation. Define 



F^v = Tr {0G M „} - -Tr {</> [D^, D v </>]} 



(11) 



Ffj. v is gauge invariant and color singlet, and such are separately the two terms 
in its definition. 

Theorem [IB] . A necessary and sufficient condition for the cancellation of 
bilinear terms A^A^ between the two terms in the right hand side of eq.(ll) is 
that 

$ = = U{x)^ dlag U\x) (12) 
with U(x) an arbitrary gauge transformation and 



( " t N Z a \ 

N — a N — a a a 



N 



N 



V 



(13) 
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For any choice of the form eq(12) F^ v obeys Bianchi identities and the iden- 
tity holds 

Flu = BpTrWAJS - d u Tr{<f> a A^} - -Tr{<^ [8^,8^]} (14) 

Ffui is gauge invariant and can be computed in the gauge in which $ a is 
diagonal. In that gauge 

F% = d^Tr{<5> a A y ) - d v Tr{^ a A^) (15) 
has an abelian form. By developing A^ iag in terms of roots A^ — a l A 1 ^ 

a 1 = diag(0, 0, ... 1,-1, ... 0) (16) 

with 

/Y;o'"!", = 5 ij 

F* v = dp A® - d u A; (17) 

The gauge transformation U(x) which brings to the unitary gauge is called 
abelian projection. 

A magnetic current can be defined as 



J a — 8 F a * 



(18) 



This current is identically zero due to Bianchi identities, but can be non zero 
in compact formulations like lattice regularization , in which Dirac strings are 
invisible. In any case J° is identically conserved 



d M J£ = 



(19) 
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and defines a magnetic U(l) conserved charge. This magnetic U(l) symme- 
try can either be Higgs broken, and then the system is a magnetic ( dual) 
superconductor, or it can be realized a la Wigner,and then magnetic charge 
is superselected. For any choice of the field $ a in eq (12 ) a magnetic U(l) 
symmetry is defined. 

To detect dual superconductivity the vev of a magnetically charged operator 
can be used as an order parameter. Such an operator has been constructed 
^H],and is magnetically charged and U(l) gauge invariant ■ The continuum 
version of the construction goes as follows. Define 

V a (x,t) = e *fd 3 yTr(r£{y,t))b ± (S-y) (20) 

where 4> a is defined by eq(12) and E(x,t) is the chromoelectric field operator 
Ei = G 0i and 

V61 = , V A 6j_ = — + Dirac string (21) 
9 r A 

[i a is gauge invariant by construction if <f> a transforms in the adjoint repre- 
sentation. In the abelian projected gauge, where $° = § a diag it assumes the 
form 

/i a (x,i) =exp J d 3 yE a ± (y,t)b ± (x-y)^ (22) 

where Ej_ is the component of the electric field along the residual U(l) 
direction as defined by eq.(16),(17), and only the transverse part survives in 
the convolution with b± . In any quantization procedure E°^_ is the conjugate 
momentum to A a ± so that fi a is nothing but the translation operator of Aj_ and 

fx a (x,t)\A a ± (y,t)) = \l a ± (y,t)+b ± (x-g)) (23) 

[i a creates a magnetic monopole. 

In the confined phase , in which monopoles condense and the ground state 
is not an eigenstate of the magnetic charge (fj, a ) ^ can signal dual supercon- 
ductivity. In the deconfined phase (/i a ) = All this refers to a given choice of 
the abelian projection, i.e. of the gauge transfomation U(x) defining . 

To explore 20 12 lj how physics depends on the choice of the abelian projec- 
tion let us go back to eq (20). By use of the cyclic invariance of the trace \i a 
can be rewritten 




In computing the correlation functions of /x a 's a change of variables can be 
performed in the Feynman integral corresponding to a gauge transformation 
generated by U(x). If U(x) is independent of the field configuration the jaco- 
bian of the transformation is 1 , the operator /i a assumes to all effects the form 
eq.(24) so that the correlators , and in particular the one point function (/j, a ) 
are independent of U(x). If U(x) depends on the field configuration , as hap- 
pens e.g. for the max abelian gauge or for any gauge in which a specific field 
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Figure 2: An example of probability distribution of the difference of the two 
highest eigenvalues of the phase <E> of the Polyakov line e**, at the lattice sites. 
SU(3) gauge group, (5 — 6.4, lattice 16 4 , 10 3 configurations. 

dependent operator is diagonalized, then the jacobian can be different from 1 
and the correlators depend on the abelian projection. However, if the number 
density of monopoles is finite ,the gauge transformation which connects two 
abelian projections is continuous everywhere except in a finite number of points 
and preserves topology: the operator ^ a defined by eq(22) will then create a 
monopole in all abelian projections. If (/i a ) ^ it signals dual superconductivity 
in all abelian projections. 

An extensive investigation of the density of monopoles in different abelian 
projections has been performed, and indeed the number density of monopoles 
is finite. Fig 2 illustrates the method , and refers specifically to the abelian 
projection in which the Polyakov line operator is diagonal. The eigenvalues of 
that (unitary) operator have the form 



and in defining the abelian projection are ordered in decreasing order of <f>i. A 
monopole singularity in a point x implies that two eigenvalues are equal , eg <pi 
and 4>2 ■ Fig 2 shows the distribution of the difference of the first two eigenvalues 
on the lattice sites of 1000 field configurations of quenched SU(3) on a 16 4 lattice 
. In no site there is a monopole. Repeating the determination on a finer lattice 
gives similar results. As a consequence one can state that the number density 
of monopoles is finite and the dual superconductivity (or non) is an intrinsic 
property , independent of the abelian projection which defines the monopoles. 

An extensive analysis on the lattice JHj shows that the vacuum is indeed a 
dual superconductor in the confined phase, and goes to normal in the deconfined 




i = 1,2,3 



(25) 
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one. A finite size scaling analysis of the susceptibility p denned as 

P° = ±W) (26) 

gives that in the quenched case {/i a ) is strictly zero above the critical tempera- 
ture defined by the Polyakov line order parameter, is different from zero below 
it. The behaviour around T c allows to determine the critical indices , which 
are consistent with those determined by use of the Polyakov line. This is clear 
evidence that dual superconductivity is a mechanism for confinement 1.8; ■ 



4 The case of full QCD 

The order parameters (p a ) can equally well be defined in the presence of dynam- 
ical quarks (Full QCD) [23 an d have the same physical meaning of creators of 
monopoles. One can then ask if a criterion for confinement could be provided by 
(p a ) , s , i.e. by dual superconductivity (or absence of) . One should prove that 
in the dual superconducting phase no colored asymptotic states exist, which is 
of course non trivial. However this is not much different from the situation in 
quenched theory with the Polyakov line criterion , as discussed in sect 2. It 
has in fact been checked [221 that QCD vacuum is a dual superconductor in 
the phase below the critical line of figl, ((/i a ) ^ 0), and is normal in the region 
above it ((p a ) — 0). The finite size scaling analysis in this case goes as follows. 
For dimensional reasons the order parameter (/i a ) has the form 

( M °) = $(|,^,mL^) (27) 

where £ is the correlation length, a the lattice spacing, m the quark mass and 
Uh the corresponding anomalous dimension . Near the critical line £ goes large 
compared to a and the dependence on a/£ can be neglected. (Scaling) The 
problem has two scales. If yh is known one can choose different values of the 
mass and of the spacial size N s such that mN s is constant, and then 

N 

<M°> T ^/(y) (28) 

By use of eq(9) the variable iV s £ can be traded with tN^ v and the scaling law 
follows 

p a = N^f(rN^) , r = (1 - ^) (29) 

whence v can be extracted and the order of the transition can be determined. 
The result is v — .33 compatible with a first order transition. 

A cross check is obtained by studying the scaling of the maximum of the 
specific heat, which for the same choices of m and N s should scale as 

C v - C° ~ Ns a/u (30) 
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If the critical indices determined through (fi a ) coincide with those resulting 
from the analysis of the specific heat, this would be additional evidence for dual 
superconductivity as a mechanism of confinement, implying that (fi a ) can be 
the order parameter. 

The situation is described in ref[2H| ,and is presently at the stage of indication 
that this is indeed the case. Numerical work is on the way which will definitely 
clarify the problem. 

In conclusion Confinement is a fundamental but difficult problem. Some 
understanding has been reached on the symmetry patterns involved. Lattice is 
a unique tool to address the problem. 

Thanks are due to my collaborators J.M. Carmona, L.Del Debbio, M. D'Elia, 
B. Lucini, G. Paffuti, C. Pica for discussions . This work is partially supported 
by MIUR Progetto Teoria delle Interazioni Fondamentali. 
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